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ABSTRACT 

We present observations of NGC 839 made with the Wide Field Spectrograph on the ANU 2.3m telescope. 
Our data cover a region 25" x 60" at a spatial resolution of ^ 1"5 The long axis of the field is aligned with 
the superwind we have discovered in this starburst galaxy. The data cover the range of 3700-7000A, with a 
spectral resolution R ^ 7000 in the red, and R ~ 3000 in the blue. We find that the stellar component of 
the galaxy is strongly dominated by a fast rotating intermediate-age (^ 400Myr) A-type stellar population, 
while the gas is concentrated in a bi-conical polar funnel. We have generated flux distributions, emission line 
ratio diagnostics and velocity maps in both emission and absorption components. We interpret these in the 
context of a new grid of low-velocity shock models appropriate for galactic-scale outflows. These models 
are remarkably well fit to the data, providing for the first time model diagnostics for shocks in superwinds 
and strongly suggesting that shock excitation is largely responsible for the extended LINER emission in the 
outflowing gas in NGC 839. Our work may have important implications both for extended LINER emission 
seen in other galaxies, and in the interpretation of objects with "composite" spectra. Finally, we present a 
scenario for the formation of Eh-A galaxies based upon our observations of NGC 839, and its relation to M82. 

Subject headings: galaxies: active — galaxies: evolution — galaxies: individual (NGC 839) — galaxies: ISM 
— ISM: jets and outflows — shock waves 



1. INTRODUCTION 

The physical nature of the low-ionization narrow emission- 
line regions (LINERs) in galaxies has been an a rea of in- 
tense study and deb ate since their d e finition by iHeckmanI 
(IT98(ih. Work bv iHo et alJ (fT997l) . iKewlev et alJ (120061) 
and Cid Fernandes et al.l (12010 ') inter alia suggests that this 
class of object represents the most common type of non- 
starbursting galaxy. LINERs, as their name suggests, are de- 
fined by their relatively strong low-ionization state lines (e.g., 
[O I], [N II]and [S II]). They are found both in the nuclei 
of giant elliptical galaxies, or, as in the case of NGC 839, in 
lower mass gas-rich star-forming systems. 

The excitation mechanisms proposed to explain the LINER 
phenomenon have included excitation by weak or low- 



luminosity active galactic nuclei (AGNs: iFerland & Netzei 
(IT983|) : iHalnern & Steiiie^ (H^SS^; 'Ho' ('1999h: IKewlev etal. 
(120061) : iHd dZOO^TlEracleous et al. (201^^ young 
and/or po st-asymptotic giant branch (AGB) stars 
(iTerlevich & Me lnick 19851 Trinchieri& di Serego Alighieril 



1991; Shields |1992t iBinette et all I1994T Stasihska et alj 
2008; .Sar zietaP 120101). shocks (.Heckman ,198C. 
Dopita & S utherland ll996riKewlev et aLllgOO la: LiparisLal 



2004 and cooling flows (iVoit & DonahuellT997l) . Indeed it is 



likely that in many galaxies, and in different regions of the 
same galaxy, more than one of these mechanisms may be at 
work. Recent stud ie s of mostly pas s ive, e ar ly-type galaxies 
dSarzi et all 120 lOt lAnnibah et all I20IOI: lEracleous et al] 
1201 eft present evidence both for and against weak AGNs, 
post- AGB (PAGB) stars and shocks in various systems. 
Although these authors find some evidence for excitation due 
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to PAGB stars or weak AGNs, neither mechanism provides a 
wholly satisfactory answer in all cases. The most compelling 
evidence in support of the PAGB mechanism is the finding 
that (in elliptical galaxies) the emission-line flux correlates 
very well with the ho st galaxy stellar l umino sity within the 
emission-line region (^Macch etto et alj II996I) . and that the 
emission-lin e flux d istribution follows that of the stellar 
continuum (Sarzi et al. 2006, 2010"). However. 'Binett e et alj 
(1994) and Cid Fernandes et al. (2009) demonstrated fliat 
photoionization by PAGB stars can only explain LINERs with 
relatively weak emission lines. The Ha equivalent width 
produced by an 8 Gyr old and a 13 Gyr passively-evolving 
stellar population amounts to 0.6 A and 1.7 A, respectively. 
Stronger emission lines require far more ionizing photons 
than old stars can provide. 

However, studies of more active galaxies and merg- 
ers, notably those with winds, find strong evidence 
for the dominance of emission by shock ex citation 
jVeilleux & Rupke 2002; Veilleux et al. 2003; Lipa ri et al. 
2004t IMonreal-Ibero et alj 1201 ft ISharp & Bland -Hawthorr 
20101). The line ratios observed in the ionized gas are typi- 
cally consistent with LINER-like emission, although this class 
of emission can extend far from the nucleus and disk of the 
host galaxy. In this case, the gas may lie far from the nucleus 
and is not ionized directly by nuclear activity. In more dis- 
tant galaxies, the entrance aperture of the spectrograph may 
capture both the circum-nuclear and the extra-nuclear compo- 
nents, making the classification still more difficult. 

In order to separate the various spatial components of 
the ionized plasma, and to disentangle the potentially com- 
plex kinematics and diverse array of power sources in galax- 
ies, we need to employ integral field units (IFUs). These 
are increasingly providing insight into complex galactic sys- 
tems. Several wider-field IFUs have been commissioned 
and u sed, notably the SAURON instrument dde Zeeuw et alj 
120021) . which has been extensively used to survey early-type 
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galaxi es. The SINS surv ey employed the SINFONI instru- 
ment ( IGenzeletal J 120081) to exam ine high-z systems. The 
VIMOS IFU dArribasetal] 120081) has been used for a sur- 
vey of Ultraluminous In frared galaxies, and th e PINGS survey 
of nearby disk galaxies ("Rosales-Ortega et al. 2010). Finally, 
the SPIRAL instrument (Sharp & Bland-Hawthorn 2010) has 
been employed to look at galactic winds. In this paper, we will 
describe data obtained of the superwind galaxy NGC 839 us- 
ing the Wide Field S pectrograph (WiFeS). Th is instrument is 
described in detail in lDopita et aLl (l2007i EoTOl) . 

NGC 839 is a nearly edge-on disk galaxy (/ = 67°, PA.= 
350°) in a Hickson Compact Group (HCG 16): a dense, in- 
teracting group of seven spi r als, six of which are actively star 
forming (JRubin et al.lll"99lt iRibeiro et al.l[T996l) . NGC 839 
may well be a recent addition to the group, more kinemati- 
cally distu rbed and showing signs of an ongoing or recent in- 
teraction (Rubin et al. 1991-. Mendes de Oliveira et al.' fl998l: 
Ide Carvalho & Coziol 1999). Belsole et al. (2003) have de- 
tected diffuse, hot X-ray emitting gas encompassing the whole 
group, evidence that it forms a bound group even though it has 
no bright early-type galaxies. 

Multiwavelength observations of NGC 839 indicate 
both ongoing star formation and LINER-like line ratios. 
iRibeiro et al.l (11996) used the [N II]/Ha ratio to classify 
NGC 839 as a starburst, though their measured values ac- 
tual ly fall in the comp osite starburst/ AGNs region when 
the 'Kewlev et al.' ('2006') scheme is applied. In addition, 
lie Carvalho & Coziol (1999) have employed the observed 
[N II]/Haand [S II]/Ha ratios to classify NGC 839 as a 
LINER + Seyfert double nucleus. X-ray spectra of the galaxy 
indicate the presence of an active starburst coupled with 
slightly higher than solar metallicities and may indicate the 
prese nc e of a low luminosity, obscur ed AGN dTurner et al.l 
120011) . iGonzalez-Martin et all (l2006h however find no evi- 
dence for an unresolved point source in hard X-ray images 
and classify NGC 839 as a starburst galaxy using X-ray mor- 
phology. 

At IR-wavelengths, NGC 839 is clearly classified as a 
Luminous Infrared Galaxy (LIRG) with Lju ^ IO^^Lq 
(lArmus et a l. 2009). This luminosity derives from a nu- 
clear starburst, which is probably driving an M 82-like galac- 
tic wind, evidenced by the very extended polar funnels of 
Ha emission seen by the Surve y for Ionization in N eutral Gas 
Galaxies (SINGG) Ha images (iMeurer et al.l(l2006<') : see Fig- 
ure [T]!. 

The observations presented in this paper are the first results 
from the WiFeS GOALS survey. The WiFeS GOALS survey 
is a sub-sample drawn from the larger Great Observatory All- 
Sky LIRG Survey (GOALS) sample. We discuss the observa- 
tions and data reduction in Section|2] In Section[3] we present 
the results of our observations and analysis and present spec- 
tra, rotation curves, and line ratio maps. In Section IH we 
compare new photoionization models and low-velocity shock 
models to the observations, and compare the derived proper- 
ties of NGC 839 to the prototypical superwind galaxy M82. 
We conclude that the balance of evidence is in favor of the su- 
perwind of NGC 839 being shock-excited. Finally, in Section 
|5] we discuss the relationship of NGC 839 to "composite " 
(or transition) galaxies and suggest extended "composite " or 
LINER emission may be fueled by shock emission. We also 
suggest that NGC 839 is one possible "missing link" to Eh-A 
galaxies and suggest an evolutionary scenario to explain how 
this class of objects might have arisen. 

Throughout this paper, we adopt the cosmological pa- 




FlG. 1. — Co mposite R+Ha (green + red) image of NGC 839 
from lMeurer et all (2006 ) and tlie archival Hubble Space Telescope Advanced 
Camera for Surveys {HST ACS) composite B+R+I blow-up of the disk re- 
gion. The two WiFeS pointings are indicated by the green rectangles overlaid 
on the SINGG image. 

rameters i?o=70.5 km s"^Mpc~\ rj v=0.73, and Om =0.27, 
based on the 5 year WMAP results dHinshaw et alj 12009). 
With these parameters, and taking into effect both the Virgo 
and Shapley infall, the adopted distance to the HCG 16 is 
55 ± 1 Mpc (NASA Extragalactic Database, NED). This 
gives a spatial scale for NGC 839 of 260 pc arcsec^^. 



2. OBSERVATIONS, DATA REDUCTION AND ANALYSIS 
2.1. Observations 

Our data of NGC 839 were taken using the WiFeS at the 
Mount Stromlo and Siding Spring Observatory 2.3 m tele- 
scope. WiFeS is a new, dual beam, image -slicing IFU de- 
scribed in detail by'Dopita et al.ld2007ll20 10h. The instrument 
is mounted in a stationary position at a Nasmyth focus to pro- 
vide excellent stability. The IFU field consists of 25 x 1" wide 
slitlets, each of which is 38" long. The spatial pixel is 0".5 
along the slitlet axis and 1".G in the spectral direction. In the 
blue, the spectra cover the range of 3700-5700 A, at a spectral 
resolution of i? ~ 3000 (« 100 km s~^). In the red, we cover 
5700-7000 A at a resolution of i? - 7000 (« 45 km s^^). 
Thus, the data have sufficient spectral resolution in the red to 
allow detailed dynamical studies, while the total wavelength 
coverage is sufficient to support excitation and chemical abun- 
dance analyses. 

We observed NGC 839 on 2009 August 18, about 4 months 
after the instrument had been commissioned. Two pointings 
were made at a position angle (PA.) of 350° (see Figure[T]), so 
that the slitlets were aligned with the minor axis of NGC 839. 
This allows us to best probe the previously detected extended 
Ha emission which occurs in two funnels both abo ve and 
below the disk (.Mendes de Oliveira et al.i il998: M eurer et alj 



120061). Each pointing consisted of three 800s exposures giv- 
ing a total integration time of 2400s for the extended emission, 
and 4800s in the overlapping region centered on the disk. The 
seeing was ^ 1.2" during the course of our exposures. Sky 
subtraction was carried out using a sky spectrum drawn from 
relatively clean sky regions at the ends of each data cube. 

2.2. Data Reduction 

The data were reduced an d flux calibrated usin g the WiFeS 
pipeline, briefly described in'Dopit aet al.l (|2010'), which uses 
IRAF routines adapted primarily from the Gemini NIFS data 
reduction package. 

The bias subtraction is somewhat complicated by the use of 
quad-readout to decrease chip read-out time as well as due to a 
slight slope and instability in the bias across each region of the 
chip. Bias frames are taken immediately before and after each 
set of observations and a two-dimensional fit of the surface is 
subtracted from the temporally nearest object data in order 
to avoid adding additional noise to the data. Any resulting 
residual is accounted for with a fit to unexposed regions of 
the detector. 

Quartz lamp flats are used to account for the response curve 
of the chip and twilight sky flats are used to correct for illu- 
mination variation along each of the slitlets. Spatial calibra- 
tion is carried out by placing a thin wire in the filter wheel 
and illuminating the slitlet array with a continuum lamp. This 
procedure defines the center of each slitlet. The individual 
spectra have no spatial distortion because the camera corrects 
the small amount of distortion introduced by the spectrograph. 
Thus only low-order mapping of the slitlets is required. 

Wavelength calibration is performed using a CuAr arc lamp 
to provide sufficient lines in both the blue and red arms of the 
camera. Arc lamp data were taken in between sets of ob- 
ject exposures. Each of the 25 slitlets is then rectified by the 
pipeline into a full data cube (one for each arm) sampled on 
a common wavelength scale. The resulting data cubes from 
each exposure were flux and telluric calibrated using observa- 
tions of the white dwarf Feige 110 and corrected for the effect 
of atmospheric dispersion. 

Cosmic -ray removal was performed with the "dcr" routine 
([Pych 2004). The final reduced and flux-calibrated data cubes 
were binned by 2 pixels along the slit in order to increase 
signal-to-noise ratio and to produce square spatial elements 
1" X 1", corresponding to ^ 260 pc at the distance of the 
galaxy. 

We analyzed each spectrum using an automated fitting rou- 
tine written in IDL, UHSPECFIT , which is based on the 
code created bv IZahid et al.l (120 lOJ) and is also employed by 
iRupke et al.l (120101) . Our routine fits and subtracts a stellar 
continuum f rom each spectrum usin g population synthesis 
models from lGonzalez Delgado et aTl ^2005) and an IDL rou- 
tine which fits a linear combination of stellar templates to a 
alaxy spectrum using the method of Moust akas & KennicuttI 
2006). Although our population synthesis models include 
stellar absorption features, it was necessary to fit excess ab- 
sorption in H/3or Ha in some spectra to ensure a reliable 
emission line flux estimates. We believe this is due to the 
inadequacy of the stellar models. Lines in the resulting emis- 
sion spectra are fit using a one or two-component Gaussian, 
depending on the goodness of fit determined by the routine. 
All of the emission lines are fit simultaneously using the 
same Gaussian component or components. Resulting fits were 
given a cursory inspection by eye to ensure the routine did 
not fail. Both continuum and emission lines were fit using the 




Fig. 2. — Example of the line and continuum fitting procedure described in 
the text applied to two spatial pixels extracted from the image cube. Note the 
clear line splitting seen in the second example. Each box has been scaled to a 
total velocity range of 3000 km s^^ in the x-direction and from an intensity 
of zero up to the peak intensity in the {/-direction. The raw data are in black, 
the continuum fit is in green, and continuum+emission fit is in red. The lower 
blue curve gives the residuals to the fit. 

MPFIT package, which performs a least- squares analysis u s- 
ing the Levenberg-Marquardt algorithm (Markwardt 2009). 

The fitting routine provides a redshift, flux, and width for 
each Gaussian component in each spatial pixel. The blue and 
red spectra are fit simultaneously, while accounting for slight 
variations in wavelength zero point and instrumental resolu- 
tion between the two arms as parameters in the fitting routine. 
The H/3 , [O III] 5007 A, [O I] 6300A, [N II] 6583A, Ha , 
and [S II] 67 1 7,673 1 A emission lines were fit reliably across 
the disk and in the extended emission. Where there are two 
components, we consistently apply the "principal" label to the 
component with the strongest Ha emission line flux. 

We apply a signal-to-noise cutoff of 5a for accurate mea- 
surement of emission line fluxes. This sets an effective flux 
sensitivity of ^ 2.5 x lO^^^erg s^^ cm^^ in [O I], [N II], 
[S II] and Ha and - 1.5 x IQ-^^erg s~i cm'^ in [O III] 
and H/3 . In all, some 300 spatial pixels (spaxels) lie above 
the > 5ct cutoff in the strongest line ([N II] or Ha ) out of a 
total of ^ 1500 possible spatial elements in the final binned 
and mosaicked cube. 

An example of one of our fits is shown in Figure 2. In this 
example, the line splitting is clearly evident, as well as the un- 
derlying stellar continuum. In this example, all of the emis- 
sion lines shown in Region 1 are above the signal-to-noise cut- 
off mentioned above, while in Region 2 only Ha and [N II] 
are above the cutoff, the rest of the emission line fluxes are 
thrown out for the purposes of our analysis. The cutoff is ap- 
plied to each component in order to avoid contamination from 
a non-detected second component in weaker lines — in Figure 
2, the [O I] line in Region 1 is an example of a line where 
the primary component is above the cutoff while the second 
component is not. 

3. RESULTS 

3.1. Overall Structure 

The archival image of NGC 839 obtained with the HST 
(proposal 10787, PI: Charlton) shows a galaxy which has a 
remarkable morphological resemblance to the famous super- 
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Fig. 3. — Blue WiFeS spectrum extracted from an aperture 3". in diameter 
located in the eastern portion of the disk. The spectrum is presented in rest- 
frame coordinates. The dominant E+A-type spectrum — which is prominent 
throughout the stellar disk — clearly shows that the stellar luminosity in the 
disk of NGC 839 is dominated by an intermediate-age (~ 500 Myr) post- 
starburst population. 

wind galaxy M82 (Figures [T] and \T3[ . It displays a highly 
inclined irregular disk of stars crossed by aligned filamen- 
tary dust lanes and with a nuclear starburst region covering 
~ 1.5 kpc. Some portions of this nuclear starburst are heav- 
ily obscured by the dust lanes, particularly toward the east, 
giving the impression in low-resolution images that there is 
more than one nucleus. The images generated from the con- 
tinuum channels of our data cube show a similar structure, 
albeit of lower resolution. Our data show a nearly edge on 
disk galaxy with a similar inclinatio n and P.A. to previous ob- 
servations (1= 67°, PA= 350°; cf. iMendes de Oliveira et"al] 
[l998). Images generated in the various emission lines show 
similar structure to previous Fabry-Perot and narrowband fil- 
ter observations of the Ha gas. The line-emitting gas extends 
somewhat asymmetrically above and below the disk by sev- 
eral kpc as traced by the red Ha emission in Figure 1 . 

The HST image reveals the details of the starburst nucleus. 
The central star-forming concentration is very bright and cov- 
ers some 230 pc. It is surrounded by what appears to be a 
disk of lower-luminosity clusters, with the outer edge delin- 
eated by a complex of dust lanes. Most of the luminosity in 
the main starburst nucleus is concentrated in a region^ 850pc 
in diameter 

Our WiFeS observations show a strong stellar continuum in 
the main stellar disk dominated by an A star-like spectrum. 
An example of this stellar spectrum taken from outside the 
central line-emitting region is shown in Figure [5] This spec- 
trum, taken by itself, would place NGC 839 in the class of 
E-i-A (elliptical plus A-star) galaxies which have undergone a 
starburst 100 - 1000 Myr ago. However, the spectrum of the 
core region of the galaxy is clearly that of an active starburst, 
although a weaker underlying A-type absorption is still seen, 
suggesting that the galaxy is currently undergoing a transi- 
tion from a starburst to an E-i-A galaxy. The stellar contin- 
uum displays stronger attenuation in the eastern portion cor- 
responding to an increase in the dust extinction accompanied 
by a slight increase in electron density n^, as traced by the 
[SII]6717/6730A ratio. Elsewhere in the disk the [S II] ratio 
provides electron densities on the order of rig ~ 100 cm"'^, 
decreasing to the low density limit (rig ^ 50 cm"'^) away 
from the disk. 
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Fig. 4. — Emission line gas rotation map of NGC 839 superimposed on 
the HST ACS F814W archival image. Where the Ha line profile is fit by 
multiple components, we have used the flux-weighted mean velocity of all 
components to construct this curve. The rotation curve has not been extended 
to the polar regions because these are dominated by the outflow motions in 
the superwind. The maximum velocity of rotation is 100 kms~^ and we infer 
that the enclosed mass in the central ~ 1.5kpc is 1.9 it 0.2 X 10^ Mq. 
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Fig. 5. — Rotation curve generated from the ionized gas velocity (red 
points) and stellar velocity (blue points). Values plotted are extracted from 
a 2" wide region along the major axis (P.A.= 350°). Average error bars for 
the gas and stellar components are given. Gas velocities are an average of the 
principal component and stellar velocities are measured from binned spec- 
tra. The stellar rotation curve traces the gas rotation within the eiTors, though 
there may be an indication of flattening in the eastern part of the disk beyond 
the ionized gas. Each spatial pixel corresponds to l"or 260 parsecs at the 
distances of NGC 839. 

3.2. Rotation Curve and Mass 

The majority of the emission line gas can be fit by two 
Gaussians. Along the major axis of the galaxy, we find a gas 
rotation curve that is more or less well behaved in the inner re- 
gions, but which appears to turn down again further out. Fig- 
ure |4] shows the rotation of the ionized gas, superimposed on 
the HST archive image, the apparent maximum in the gas rota- 
tion is coincident with the dust lanes. Figure |5] shows rotation 
curves of the ionized gas traced by emission line velocities 
and stars traced by Balmer absorption, extracted from a re- 
gion 2"wide along the major axis. The stellar data have been 
averaged from the velocities measured for each spaxel, the 
stellar data are binned, and the principal component is plotted 
here. The stellar rotation is more difficult to measure in the 
starburst region due to weak stellar absorption features, but 
appears to trace the gas rotation. Beyond the central starburst 
the stellar rotation may flatten, though it is difficult to discern 
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Fig. 6. — Here, we compare (a) the nebular and stellar rotation curve at H/9 
with (b) the He I5876A emission line and much more complex interstellar 
Na D absorption lines. Each panel integrates a 3"diameter aperture, and are 
separated by 3"each along the major axis of NGC 839. The short bars on 
each panel give the local systemic velocity of the emission line gas for each 
spectral feature. Again, the predominance of the A-type stellar absorption is 
prevalent in the outer regions of the galaxy. Note that the stellar component 
has a faster rotation velocity than the ionized gas in the outer regions (E3 and 
W3). The interstellar Na I lines display a much more complex structure with 
strong line splitting. One component follows the rotation of the ionized gas, 
while the other is blueshifted. We infer that this arises in a neutral outflow at 
the base of the superwind with a projected velocity of ~ 130km s~ ^ . 

this given our errors. The maximum in the gas rotation at a 
radius of 750 pc corresponds to the boundary of the central 
starburst region and the enclosed mass Mxsin^(z) estimated 
to lie within this radius isl.9±0.2xl0^ M© . The peak of 
the rotation curv e, at just over 100 km s~^, is similar to that 
found by others (Ru bin et al.lll991l; iMendes de Oliveira et alj 
[l998). The maximum rotational velocity is very similar t o 
that found in M 82 by IShopbell & Bland-Hawthoriil (I1998I) . 
However, the region enclosed is 50% larger and thus the en- 
closed mass correspondingly larger. We will make a detailed 
comparison of NGC 839 with M 82 in Section |43] 

The majority of the emission line gas can be fit by two 
Gaussians. In the plane of the galaxy, as remarked above, 
the gas shows an organized rotation curve. This is also true 
of the stellar component, shown in Figure |6] Here, we have 
extracted the spectrum in the region around H/3 using sev- 
eral apertures along the major axis of the galaxy each 3"in 
diameter This figure also clearly shows the dominance of 
the A-type stellar absorption component in the outer parts of 
the galaxy. The gas appears to have a maximum in its veloc- 
ity coincident with the dust lanes apparent in the HST image 
(Figures|4]and|5]l. 




Fig. 7. — Normalized kinetic energy content -Fhq""'^ of the ionized gas 
along the line of sight, scaled to the peak value. The coordinate scale is in 
arcsec with the nucleus at ~ [11, 12.5], and the analysis has been carried 
out on individual WiFeS pixels (see Figure [T). The four boxes refer to the 
energy content of the following components: (a) the energy associated with 
the velocity dispersion of the principal line component, (b) the energy associ- 
ated with the velocity dispersion of the secondary component of the line, (c) 
the energy content associated with the line splitting, and (d) the sum of these 
three components. Note how well the outflow axis is defined by component 
(c), as would be expected for an expanding bicone of gas. 



3.3. Outflow 

The extended Ha morphology revealed by the SINGG 
images and by maps generated from our data is typical of 
that seen in gala xies w it h bico nical galactic scale superwind_s 
(e.g., Heckman etalJ 119871 1990; Veilleux et al. 199 
Lehnert & Heckman 1996; Shopbell & Bland-Hawthori 
1998; Veilleux & Rupke 2002). The emission extends 
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enussion extends a 
few kpc above and below the disk along the minor axis in a 
conical shape, though due to the distance of NGC 839 and 
the available ground-based data it is difficult to determine 
detailed characteristics of the bicone. 

Our spectra also reveal kinematic evidence for a large-scale 
outflow or superwind extending both above and below the 
disk. We see extensive line splitting with typical Av ^ 250 
km s^^ along the minor axis up to ^ 2 kpc above and below 
the plane. The maximum velocity splitting seen at any point 
in the flow is Av ^ 305 km s^^. 

As evidenced by Figure |6] there is also evidence for a 
blueshifted Na D component, which traces the neutral compo- 
nent of the superwind. The largest separation in this NGC 839 
is inclined such that we observe the blueshifted component 
below the disk, along lines of sight toward the disk itself. As 
the stellar continuum drops rapidly below the disk, we are un- 
able to trace the extent of the Na D much beyond the disk and 
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so focus mainly on the emission hne data. 

The velocity dispersion of the emission line gas is gen- 
erally lowest in the disk, on the order of 50 km s~^, and 
increases above and below the disk to 100-200 km s^^, 
consi stent with the kinematic signatures of a galactic wind 
(e.g.. Lehnert & Heck man 1996). Where two components are 
fit the ratio of the velocity dispersions are within a factor of 
~ 1-2.5 of one another Tracing the velocity dispersion is 
somewhat complicated due to the presence of multiple Gaus- 
sian components of varying relative flux and is limited by our 
spatial and spectral resolution, so in order to quantify these 
signatures, we have analyzed the data in terms of the line flux 
and the square of the velocity dispersion, u^, along the line 
of sight. The Ha line flux is given by Fho oc aBnuUcAR, 
where a is the effective Case-B recombination coefficient, nn 



is the hydrogen density in cm 



-3 



is the electron density in 



cm^"^, and AR is the path length through the ionized medium. 
This equation can be expressed more simply as Fhq cx J^Uc, 
where S is the ionized mass surface density. Thus to the ex- 
tent that the electron density is relatively constant, the Ha 
flux is a measure of mass surface density. Therefore, the prod- 
uct FiiaO'^ is a measure of the kinetic energy content of the 
ionized gas along the line of sight. The total energy associated 
with line splitting is calculated by substituting line width with 
the deviation of the central peak of each component from the 
flux-weighted average velocity and summing the product of 
each components 'deviation with the square of its flux. 

In Figure |2] we plot the various components of this energy 
content, namely, the energy associated with the velocity dis- 
persion of the principal and secondary components (panels 
(a) and (b)), the energy content associated with the line split- 
ting (panel (c)), and the sum of these three components (panel 
(d)). In the figure, the energy scale has been normalized to 
the largest value in the data cube, so as to facilitate compari- 
son between the components. 

The energy associated with the velocity dispersion in each 
of the components shows a marked tendency to be high within 
the two funnels of ionized gas extending almost perpendicular 
to the stellar disk. This alignment is much more evident in 
terms of the energy content associated with the line splitting 
shown in panel (c) and clearly defines an axis for the outflow 
(P.A.= +3°±3°). This alignment would be expected if the gas 
motions lie along the surface of an expanding funnel, since 
the greatest expansion velocities seen in projection would lie 
along the axis of the funnel in this geometry. 

Figure |7] clearly shows that we are observing an ionized 
superwind moving out of the galaxy. The space velocities 
associated with the outflow are certainly greater than the ~ 
130km s^^ implied by either the blueshift in the Na D -line 
absorption or the ~ 125km s^^ implied by the line splitting 
observed in the funnel. The Ha image of Figure [T] limits the 
cone opening angle to not more than 60°. Therefore, if the 
predominant gas flow is along the surface of the cone, the 
actual outflow velocity would be ^ 250km s^^. 

The outflow at P.A.= +3°±3° is neither perfectly aligned 
along the photometric minor axis of the galaxy or with 
the rotation axis. The photometric line of nodes measured 
by iMendes de Oliveira et al.l ([1998) for the stellar disk is at 
P.A. = 84°±5°. From die HST image, we estimate P.A. 
— 80°±3°. Thus, the photometric minor axis is at P.A.= 
352°±4°. From our rotation curve (see Figure |4|i, we in- 
fer that the polar direction is at P.A.= 342° ±3°. The mild 
disagreement between the photometric and kinematic hne of 
nodes had been previously noted bV iMendes de Oliveira et al.l 
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Fig. 8. — Hq surface brightness as a function of distance (r) from the nu- 
cleus of NGC 839. The surface brightness of the northern cone is plotted 
as blue symbols and of the southern cone as red symbols. We have taken a 
model in which the outflow has the form of two cones with opening angle 
30°meeting at the midplane. and covering the starburst region, so that the 
starburst is effectively 0.78 kpc from the apex of the cones. The solid line is 
not a fit, but slope r~^. The dotted lines show the shallower r~^ and steeper 
r~* slopes. 



(2003). This perturbation could be due to the ongoing merger 
activity in HCG 16. 

3.4. Surface Brightness 

The observations described in the previous section suggest 
that the geometry of the superwind is close to that of a hol- 
low bi conical outflow. Models of superwinds bvl Cooper et alj 
(120091) show that the optically emitting filaments are entrained 
and accelerated by a faster and hotter outflow from the nu- 
clear starburst region. In this model, the pressure of the opti- 
cally emitting material is determined by the ram and thermal 
pressure of the surrounding hot medium. The filaments are 
ionized either by the escape of UV photons from the central 
starburst, or, if the outflowing medium is optically thick to the 
escape of Lyman continuum photons, then by cloud shocks 
which propagate into the dense filaments as they are acceler- 
ated in the outflow. 

A study of the radial variation of surface brightness within 
the cone can help us to distinguish between these possibil- 
ities. Let us consider two cases, either (1) that the outflow 
has constant velocity or (2) that the outflow velocity increases 
proportional to the radial distance. Recall first the equation 
derived in the previous section; F^a oc Sno, where E is the 
ionized mass surface density and ric is the electron density. 
Since the gas pressure is given by P = nkT and the tem- 
perature of the ionized plasma is of order lO^K, we can write 
the Ha flux as Fhq oc SP. In a photoionized entrained out- 
flow, the ionized mass is conserved so for constant velocity 
of expansion, S oc r^^, where r is the radial distance from 
the cone apex. The ram pressure in the hot flow decreases 
as P oc r^^. Therefore, in a constant- velocity photoionized 
entrained outflow we expect Pho oc r^^. In case (2) (an ac- 
celerating outflow), Phq oc r~*. 

In a single shocked cloud, the luminosity at Ha is propor- 
tional to the energy flux across the cloud shock; S'c,Ha oc 
AcUcV^, where Ac is the cross-sectional area of the cloud. 
Tic is the density of the clouds (or filaments), and Vc is the 
cloud shock velocity, given in terms of the external pressure, 
P, by P = ricVc- Therefore, we can rewrite the Ha luminos- 
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ity produced by a single cloud/ filament as Sc.Ha ex AcPvc- 
In a conical outflow, A oc r^, P oc r^^, and Vc is constant. 
Therefore, the flux produced in a single cloud is independent 
of its distance from the center of the flow. However, the to- 
tal number of clouds is conserved and their space density de- 
creases as r~^ for constant velocity of expansion. In this 
case, then Fhq oc r~^. In case (2) of an acc elerating out- 
flow, f H» oc r~^. In hydrodynamic models by ICooper et al.l 
(l2009h . the outflowing filaments are shredded and destroyed 
as they are accelerated, so in this case the radial dependence 
of the surface brightness could be even steeper than r^^. 

In Figure [8] we show the variation of surface brightness in 
Ha along the biconical outflow. We have presented the mea- 
surements using a simple bicone with an opening angle of 30°, 
which i s roughly appr oximated using our data and the Ha im- 
ages of iMeurer et al. ' (2006). We assume that the cone origi- 
nates in the central starburst, which has a radius of ^ 0.26 kpc 
in the disk plane corresponding to the size of the central star- 
forming concentration seen in the HST images. In this case, 
the starburst lies 0.78 kpc from the origin of the coordinates 
for each half of the bicone. Figure [8] shows that the slope of 
the surface brightness curve is close to r^^, which would tend 
to support the case for a photoionized outflow within a narrow 
cone in the polar direction, though this does not rule out the 
presence of other ionizing sources. 

3.5. Emission Line Ratios 

We examine line ratio diagnostic diagrams which use 
[N 11]/ Ha,[S II]/Haor [O I]/Ha ratios against the 
[O III1/H/3 ra ti o. T hese diagrams were first pr oposed by 
[Baldwin et aP (119811) and IVeiUeux & Osterbrockl (I1987h to 
determine the likely ionizing source of emission line gas in 
galaxies. The ionizing source classification on these dia grams 
was subsequently updat ed by .Kewlev et al.l (1200 Ibi 1200 6) and 
iKauffmann et al.l (l2003h . These line ratios have the advantage 
that they can still be used even when reddening corrections 
are quite uncertain. 



In the case of NGC 839, we find a very systematic behavior 
in the spatial distribution in the line ratios involving the low 
ionization species as shown in Figure |9] Figure |9] shows that 
the [N II]/Ha,[S 11]/ Ha, and [O I]/Ha ratios are systemati- 
cally lower near the nucleus, somewhat low along the axis of 
the bicone, and high along the edges of the bicone. We will 
compare the observations with detailed models on the optical 
diagnostic diagrams in the following section. 

4. ANALYSIS 

Our results of NGC 839 present unequivocal evidence that 
within the last Gyr it has undergone a long-lived and massive 
starburst. Not only is the disk dominated by an E-i-A spec- 
trum originating from this starburst, but the starburst is clearly 
ongoing and is currently actively powering a superwind evi- 
denced by large bulk motions and line splitting. 

The basic physics of superwinds po wered by starbursts is 
rather simple, and has been described by Chevalier & Oegerlg 
(119791) . To power a wind, the energy injection per unit area 
into a fraction, /, of the disk gas mu st exceed its binding en- 
ergy. As lChevalier & Oegerld (119790 have shown, this means 
that the gas escaping into the halo has an initial temperature 
greater than some cr i tical v alue related to the binding energy. 
Indeed lCooper et aP (120091) have re-iterated how similar a su- 
perwind is to a simple stellar wind in an inhomogeneous in- 
terstellar medium (ISM). 

The physical parameters which control the production of 
a gala ctic wind can easily be derived, as shown by Dopitg 
(120081) . and we repeat his argument verbatim here. The input 
energy E^a per unit area is proportional to the surface rate of 
star formation S, integrated over the lifetime of the starburst, 

■^sb; 



Ein = aS^rsB. 



The escape energy of the gas is 



Ep, 



0.5/EgasWesci 



(1) 



(2) 
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Fig. 10. — BPT diagrams for NGC 839. Filled circles are the line ratios measured in spectra from the individual spatial elements, with our nuclear values 
marked by a red cross (x). Only data with a minimum of 5a in the relevant diagnostic lines for each diagram are plotted, 91 for [N II]/Ha , 87 for [S II]/Ha 
and 71 for [O I]/Ha . The region of the galaxy from which these data are taken corresponds roughly to the [O I]/Hq map in Figure |9] We have over-plotted 
the theoretica l values derived from the photoionization models described in the te xt. The plots also sh ow the earlier maximum starburst envelope as defined by 
IKewlev et alj |200 1 a) and the division between Seyferts and LINERs as defined in lKewlev et alj 420060 . 

where Egas is the surface density of gas in the disk. The star 
formation rate is related to the gas surface density through the 
iKennicuti (fl998.) star formation law 
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Thus, a galactic wind becomes possible when 
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const. 
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Therefore, a wind is favored when the escape velocity is low 
(such as will be the case in small galaxies), when the starburst 
continues for a significant period (which is not surprising, 
since long-lived star formation events contribute more energy 
to removal of gas), and when the gas surface density is high 
because the energy input per unit mass of ISM is higher in this 
case. Our observations clearly demonstrate that NGC 839 ful- 
fills all three of these conditions. Given that shock excitation 
of the filaments in the superwind is implied by the theoreti- 
cal models of this type of flow, and also by the observed bulk 
motions, line ratios, and radial dependence of surface bright- 
ness, we clearly need to investigate shock models in more de- 
tail. The shock velocities which are suggested by observation 
(100-250 km s~ ^) are much lower than the fast shock mod- 
els described bv lDopita & Suth erland (1996) and Allen et al.l 
(I2OO8) . Indeed, low-velocity shoc k emi ssion has been a much 
negle c ted field since the wo rk of ICoxl (11 972). Dopita (119761 
[T977h . lRavmondl(IT979h and lShull & McKee. (.1979.) . 

However, it is also possible to have a superwind which 
is predominantly photoionized by the escape of UV photons 
from the central starburst, and this possibility also needs thor- 
ough investigation. Indeed the behavior of the line ratios in 
Figure |9] might well be explicable in the context of a pho- 
toionization model of the cones, since we would expect to see 
relatively low-excitation gas toward the nucleus and along the 
bicone axis, while the higher excitation gas would be found 
close to the bounding ionization front along the walls of the 
bicone. Likewise, the slope of the radial dependence of sur- 
face brightness given in Figure [8] is closer to that expected 
from photoionization models. In the following subsections, 
we investigate both shock and photoionization models. 
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Fig. 1 1 . — Ratio of total shock luminosity to Ha luminosity as a function 
of shock velocity and hydrogen pre-ionization fraction as marked on each of 
the curves. 

4.1. Photoionization Models 

In a photoionization model of the superwind, we would 
expect to see relatively low-excitation gas toward the nu- 
cleus and along the bicone axis, while the higher excitation 
gas would be found close to the bounding ionization front 
along the walls of the bicone. This can be seen qualita- 
tively in Figure |9] However, this simple model does not 
stand up to closer scrutiny. As shown in Figure [TO] the 
vast majority of the gas falls within the "composite" region 
in the [N Ill/Ha diagram, which lies below the theoretical 
"maximum starburst line" defined in Kewlevetal. (2001aD 
and above the empirical line established in Kauffmann et al] 
(2003 ) for pu r e star formation, conforming to the scheme in 
IKewlev et al.l ( 120061) . This implies a radiation field with a 
harder ionizing flux than can be accounted for by pure star 
formation. We generated the photoionization models plotted 
in Figure [Jo] using the STARBURST99 and MAPPINGS IH 
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Fig. 12. — This is a blow-up of a region of the BPT diagrams for NGC 839 given in Figure [Tol Here, we have overplotted the theoretical values derived from 
the shock models described in the text. Models plotted are for Z/Zq ^2 with the exception of [N II]/Ha , where the dotted lines show the models for solar 
abundances. Individual points are in order of increasing preionization from bottom to top and lef t to right, from t o 1 in steps of 0.2. The [N II]/Ha plot also 
shows the maximum starburst (soUd black)and composite (dashed black) lines as defined bv Kewlev et alj <2001bl) : lKauffmann et alj <2003i) and Kewlev et alj 
)200fl) . The [S II]/Haand [O I]/Ha plots show the maximum starburst line (solid black) and the division between Seyferts and LINERs (dashed black), as 
defined in lKewlev et"al] (2MB) . 

and ranged in velocity from 60 to 200 km s^^ in steps of 
20 km s^^. Model grids were run using both Z/Zq = 1 
and 2. 

For a fully radiative shock propagating at a velocity Vs into 
a cloud of density ric, the surface luminosity L (erg cm~^) 
is equal to the flux of mechanical energy through the shock; 



codes (iLeitherer et al.lll999l:ISuflierland & Dopitalll993h with 



the abundance set of iGrevesse et alj (120 10). Our models are 
unable to produce hard enough ionizin g radiation to reach th e 
theore tical pure star formatio n lines of iKewlev et alj (1200 lal) . 
which llevesque et al.l (|2010l) propose is due to a lack of rota- 
tion in the steUar models. 

Our observed [S II]/Ha ratios and [O 1]/Ha ratios become 
more H Ilregion-like toward the nucleu s and predorn i nantl 
LINER-like in other regions. Although [Turner et al.l (120011 
suggest the possible presence of a low-luminosity AGN which 
could possib ly account for the H Il/com posite line ratios in 
the nucleus (Ide Carvalho & Coziollll999 ). it cannot account 
for the extended LINER-like emission. In fact, our line-ratio 
maps show a general decrease in line ratios closer to the nu- 
cleus corresponding to a stronger contribution from star for- 
mation, though our optical emission line ratios ca nnot rule out 
an obscured AGN and the X-ray ob servations ofl Tumer et alJ 
([2001) and Gonzalez - Martm et alJ (.2006.) provide somewhat 
ambiguous results. 

4.2. Shock Modeling 

We have computed a new grid of low-velocity shock models 
using the MAPPINGS I II code, an u pdated ver sion o f the code 
which was described in lSutherland & Dopita ( Il993h . In these 
models, both pre-ionization and shock velocity are treated as 
independent parameters. Allowing the pre-ionization to re- 
main a free parameter, rather than being determined by the 
computed upstream ionizing radiation field is a novel ap- 
proach which has the advantage of removing geometry de- 
pendence and showing the total range of line ratio diagnostic 
space which can be covered by low-velocity shock models. 

The shock velocities range from 60 to 200 km s^^ in 
steps of 20 km s^^, with preionization values in the range 
0.01 — 1.0 in steps of .1. The abundance set is taken 
from' Grevesse et aTl (120101) a nd we use standard (s olar region) 
dust depletion factors from iKimura et alj ( |2003i) . In using 
these depletion factors, we are implicitly assuming that the 
shock is not sufficiently fast to sputter the dust grains advected 
into the shocked region. The shock models have a transverse 
magnetic field consistent with equipartition of the thermal and 

magnetic field (B — 5/zG for Uc = 10cm~^; B (x Uc ), 



L 



^ticVi 



However, most radiation comes out in UV or 



IR lines and continuum, and so the optical luminosity does 
not reflect the shock luminosity. In order to be able to cal- 
culate the total shock luminosity, we have computed the ra- 
tio of the total to Ha luminosity, i/L(Ha) from our mod- 
els. This is given for the convenience of others in Figure 
nn The ratio L/L(Ha) asymptotes to ^ 80 at shock ve- 
locities above ~ 140 km s^^, and can be approximated as 
being proportional to the shock velocity for lower velocities; 
L/i(Ha) - 60t)3/(100kms-i). 

From our grid of shock models we have extracted the 
H/3/[0 III], [N II]/Ha,[S II]/Haand [O I]/Haline ratios so 
as to be able to construct the optical diagnostic diagrams con- 
sistent with those presented in 13. 51 These are given in Figure 
[T2I The line ratios derived from the shock models reproduce a 
significant portion of our observed line ratios from NGC 839. 
In a broader sense these line ratios provide provide a useful di- 
agnostic tool for understanding the ionizing source in galactic 
scale winds. Previous studies of galactic winds suggested the 
presence of shocks, with shock velocities consistent with our 
slow shock models (e.g. 'Sho pbell & Bland-Hawthoriill 19981 ; 
Sharp & Bland-Hawthorn 201(f^ 

The line ratios for the shock models fall in or very 
near the composite region in the [N II]/Ha diagram while 
they migrate into the LINER region in the [S IIJ/Haand 
[O I]/Ha plots. The shock model points plotted in the 
[O I]/Havs [O III]/H/3plot provide the largest separation 
in diagnostic diagram space for shock models of different 
velocities and preionizations. The shock model values for 
[O I]/Ha also trace the gas with line ratios moving slightly 
into the Seyfert region of the diagram. Although the sepa- 
ration in [N II]/Ha values derived from the shock models is 
not as significant as [S II]/Ha or [O I]/Ha , they do provide 
a larger spread between the models using different metallic- 
ities. The shock model values for [N II]/Ha corresponding 
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Fig. 13. — Side-by-side comparison of M82 (left) and NGC 839 (right). The image of M82 is an HS T ACS B,V ,I and Hq composite image, the image of 
NGC 839 is an ACS B,R and / image with the Ha map from SINGG overiaid (Mutchler et al. 2007; M eurer et alJ20 06). The resolution has been degraded for 
the M82 image to provide a better visual comparison with the much more distant NGC 839. NGC 839 has been rotated and scaled to approximately match the 
P.A. and spatial scale of M82 to provide a more accurate comparison. 



best to the observed values from NGC 839 are consistent with 
an abundance of Z/Z( n ~ 2, which i s sim ilar to the supersolar 
abundance quoted by iTurner et alJ (2001). This can be seen 
in Figure [12] where the solar abundance models are clearly 
separate in the [N II]/Ha diagram. 

The superwind in NGC 839 exhibits a more composite 
spectrum than the sho ck-dominated, starburst-dr iven winds 
examined in detail by iSharp & Bland-H awthorn' (2010). In 
the study by citetSharplO the diagnostic diagrams indicate a 
larger percentage of the extended emission further from the 
H Il-like region. In this sense the wind in NGC 839 is some- 
what more like M82, as the extended emission associated with 
the wind, while still dominated by shocks, lies nearer the H II- 
like region of the diagnostic diagrams. In contrast to M82, 
however, the extended wind emission in NGC 839 shows a 
somewhat harder radiation field with higher [O III]/H/i? ratios 
placing the extended emission in NGC 839 much more firmly 
in the composite region and bey ond t he reach of photoioniza- 
tion models as shown in Sectionl4n 



4.3. Comparison of NGC 839 with M82 

NGC 839 shows several similarities to the nearby, well- 
studied, prototypical wind galaxy M82. A visual comparison 
reveals similar morphologies, archival HST ACS images of 
both galaxies show an irregular disks crossed by dust lanes, 
with a strong central starburst (see Figure \T3[ . NGC 839 
shows more disk obscuration and less nuclear dust obscura- 
tion due to a lower inclination angle. Ha images of both 
systems show the characteristic biconical outflow, though the 
filamentary structure of the Ha is easier to discern in M82, as 
it is 10 times nearer than NGC 839. 

Table 1 shows several physical properties of M82 and 
NGC 839, which provide an even better comparison. The in- 
frared and Ha luminosities (which are associated with the 
dust and star formation) are quite comparable. The star- 
burst driving the outflows in both galaxies is strongly con- 
centrated around the nucleus. The peak of the gas rota- 
tion curve is at a slightly larger radius in NGC 839, though 
the peak velocities are nearly identical (^ 100 km s^^, 
Shopbell & Bland-Hawthorn 1998, and references therein). 
Origha et al., (,2004.) measure a gas metallicity in M82 con- 
sistent with the metalhcity of our shock models, in both cases 
super solar. 

iShopbell & Bland-HawthornI (II998I) provide line ratio 



TABLE 1 
Comparison of NGC 839 and M82 



Parameter 



M82 



NGC 839 



Distance (Mpc)" 
Mass(0)/_R'' (kpc) 

Lmiergss ^f 

Gas metallicity [O/H] •>' 



2.3 X 10^/0.75 

2 X lO"*^ 

2.3 X 10"''* 

-3.07 



55 

2.8 X 10^/0.5 

1.25 X lO"*! 

3.6 X 10'*"' 

-3.01 



" Distances from NED. 
'' Our estimate. 
"^ Heckman et al] 
■^ Meurer et al. 




Sanders et al. 1 2003 
Origha et al.. (2004) . 



maps and diagnostics of M82 that indicate stellar photoion- 
ization within the galactic superwind, though they note that 
the line ratios become more shock-like further from the 
disk, away from the starbursting region. While the spec- 
tra in NGC 839 show a more composite source of ioniz- 
ing radiation, the line ratios become more H Ilregion-like 
close to the nuclear starbursting region. This behavior is 
consistent with the line ratio structure seen in M82 and 
the conclusion of Shopbell & Bland-Hawthorn (1998) that 
close to the nucleus photoionization by the ongoing star- 
burst begins to d ominate. Finally, Ch a ndra observations 
of NGC 839 by iGonzalez-Martin et all (120061) show fun- 
nels of soft X-ray emission along the minor axis similar to 
those seen at much higher spatial re solution in Chandra im- 
ages of M82 ( S trickland et alj 12004 ). This soft X-ray gas 
is associated with shocks in the galactic wind and in the 
case of NGC 839 could be contributing very slightly to the 
di ffuse X-ray emissio n seen surrounding the entire group 
bv lBelsoleetaI](l2003l) . 

4.4. E+A Galaxy Connection 

Our observations of NGC 839 suggest a link between 
lower-mass starburst systems and Eh-A galaxies. Although 
NGC 839 is a fairly low-mass system, it nonetheless has a 
large-scale organ ized rotation, and an elliptical-h ke (r^/*) lu- 
minosity profile dMendes de Oliveira et alj[l998l) . The stellar 
population of the outer regions away from the superwind fun- 
nels is dominated by a classical Eh-A spectrum. These proper- 
ties are precisely what has been found by ,Pracv et al., (20091) 
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in a Gemini Multi-Object Spectrograph IFU study of a sam- 
ple of 10 nearby (z — 0.04 — 0.20) Eh-A galaxies selected 
from the Two Degree Field Galaxy Redshift Survey. In their 
sample, they suggest that the large fraction of fast rotators ar- 
gues against the hypothesis that equal mass mergers are the 
dominant progenitor to the Eh-A population. 

NGC 839 and its similarity with M82 suggest a simple sce- 
nario for the formation of Eh-A galaxies. In this, we start 
with a gas-rich low-mass galaxy in a cluster or loose group. 
This produces a strong interaction during a close passage with 
another galaxy. In the case of M82 this would have been 
M81, while in the case of NGC 839 this was most prob- 
ably NGC 838, which has a slightly higher I R luminosity 
QogL/LQ — 11.05 versus log L/L0 = 11.01; lArmus et all 
l2009i) . and which also shows a strong nuclear starburst, and 
evidence of tidal interaction and galactic winds in Ha im- 
ages. This interaction drives the disk gas toward the nucleus, 
triggering a strong nuclear starburst. When this occurs, the 
starburst acts as a so urce of viscosity in the accretion disk 
JCollin & Zahn|[T999D . since shocks generated by increasing 
mass loss and SN-powered shells reduce rotational angular 
momentum in the rotationally trailing ISM, and put more an- 
gular momentum into the leading ISM, encouraging the de- 
velopment of a galactic wind. The process therefore has the 
potential to run away, as gas is fed into an ever smaller and 
denser central region, increasing the specific star formation 
rate, and driving an ever stronger galactic wind along the lines 
of Equation (4). 

The rotational velocities give further evidence that gas is 
currentl y be ing fed into the central region. As shown above in 
Section 13^ the gaseous disk shows a rotational maximum at 
the boundary of the central starburst region and in the vicinity 
of the dust lanes seen in the HST image. This could be ex- 
plicable if the outer ionized disk is not in rotational support, 
and is still infalling into the central star-forming region. A 
ve ry similar rotation curve was observ ed in the case of M82 
bv lShopbell & Bland-Hawdiornl(fT998h . 



5. CONCLUSIONS 

In this paper, we have provided very strong evidence for a 
shock-excited superwind in NGC 839. Our new shock models 
cannot only reproduce the line ratios in NGC 839, but also the 
inferred velocities agree with what is derived from emission 
line and absorption line kinematics. These observations are 
also consistent with the e xpected wind-blown shock yeloci- 
ties observed by others (Shopbell & Bland-Hawthorn '1998'; 
IVeilleux & Rupke 2002; Shai-p & Bland-Hawthorn 2010), 
and confirm the model of shock-excited superwinds sug- 



gested by recent studies JSharp & Bland-HawthornI 120 lOt 
iMonreal-Ibero et £11120101) . These results are also in li n e with 
the picture advocated by ISharp & Bland-HawthornI (12010b 
who argue that shocks outshine direct stellar photoionization 
by the time a galactic scale wind is blown out, but that ongo- 
ing star formation contributes to the luminosity and ionization 
field of the galaxy close to the disk. 

Our shock models clearly demonstrate that galactic-wide 
emission from intermediate-velocity shocks can shift the 
global line ratios into the regime of LINERs and "composite" 
objects. Although in the case of NGC 839, the "composite" 
or LINER emission is generated by shocks in a superwind, we 
can envisage a number of other scenarios which could pro- 
duce extended LINER-like intermediate velocity shock emis- 
sion. These include direct galactic collisions, mergers in 
which one of the merging galaxies is falling into the stellar 
wind-blown region of another, and minor mergers involving a 
gas-rich system falling into the hot gaseous halo of an ellip- 
tical galaxy. None of these directly involve an AGN. These 
results highlight and reinforce the need for integral field data 
to interpret composite and LINER spectra. 

Finally, we have proposed that NGC 839 provides a link be- 
tween lower-mass starburst systems and the enigmatic E+A 
galaxies. Outside of the central starburst NGC 839 exhibits 
a pure E+A spectrum and even within the starbursting nu- 
cleus the underlying A-star absorption is still evident. We 
have proposed a possible formation scenario consistent with 
the history of NGC 839 given its environment and similar- 
ity to M82 whi ch agrees with the recent observational work 
of iPracv et al] (12009 ). Further integral field observations of 
systems similar to NGC 839 and M82 should provide further 
insight into our proposed idea. 
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